OITS 716 
April 2002 



Inclusive Distributions for Hadronic Collisions in the Valon-Recombination Model 

CS|"; Rudolph C. Hwa 1 and C. B. Yang 1 ' 2 

, 1 Institute of Theoretical Science and Department of Physics 

O ■ University of Oregon, Eugene, OR 97403-5203, USA 

2 Institute of Particle Physics, Hua-Zhong Normal University, Wuhan 430079, P. R. China 

5h ' 

Inclusive distributions of soft production of mesons in hadronic collisions are calculated in the 
valon-recombination model. The new determination of the valon distributions from hard scattering 
data makes possible a tightly interrelated treatment of pion and kaon production in the fragmenta- 
tion regions of proton, pion and kaon. Only one free parameter is used in the determination of the 
valon distribution in the kaon. No other adjustable parameter is needed to fit the ^-distributions 
of the data on inclusive cross sections, except for the normalizations since the data are at fixed pr- 
The success of the model in reproducing seven inclusive distributions suggests that there are two 
important mechanisms at work in soft production. One is that the structure of the hadron that 
fragments is highly relevant. The other is that the produced particles are formed by the recombi- 
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nation of quarks and antiquarks. These two aspects about hadrons in soft processes can be well 



described in the framework of the valon-recombination model 



I. INTRODUCTION 
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Particle production at low px in hadronic collisions has always been a challenge for theoretical models to describe, 
since, on the one hand, it is a process that is non-perturbative, while, on the other hand, meaningful modeling can 
only be done in the framework of quarks and gluons. The subject has a long history, some early reviews of which can 
be found in Refs. [Q, 0|. Since soft processes, as they are called, cannot be treated by perturbative QCD, they have 
not been given the degree of rigorous scrutiny that have been accorded hard processes. Among the models that are 
constructed for soft processes, there are basically only two types: string models that make use of fragmentation |3| [|] 
and parton models that are based on recombination |(|. This paper treats the modernization of the latter. What 
is new is that the parton distributions of the proton have recently been calculated by CTEQ in pQCD over wide 
ranges of x and Q 2 , fitting a large collection of experimental data [JfJ. From such parton distributions the structure 
of the proton in terms of the valons || can be more precisely determined || . With the new parameterization of the 
valon model now available, it is possible to revisit the problem of hadron production in soft processes and calculate 
the inclusive distributions of the produced particles without adjustable parameters. Good agreement with the low-pr 
experimental data can give definitive support to the recombination model. 

The kinematical region in which we focus our attention is the projectile fragmentation region, roughly x > 0.2. In 
the central region the structure of the projectile hadron is less important. But in the fragmentation region it is known 
as early as the mid-70's that the inclusive distribution of the produced pions is closely related to the structure function 
of the proton, as observed by Ochs pTOj ] ■ The recombination model || is a realization of that observation, and the 
valon model is a self-consistent formulation of the unification of hadron structure and recombination probability || . 
In short, our view is that, while the fragmentation of strings may be suitable for the central region, the recombination 
of partons is more relevant for the fragmentation region. 

It has been regarded as a striking confirmation of the dual parton model (DPM) [p| t o reproduce the charge 
distribution of the produced particles in 7r + p collisions; the forward-backward asymmetry ]ll| is presented as evidence 
for the two-chain diquark fragmentation of the proton. While the qualitative agreement with data gives support to 
the two-chain mechanism of DPM as opposed to the single-chain mechanism of the LUND model, the data (not cited 
in ||) were old and inaccurate and the theoretical calculations were crude. The data that we shall compare with 
are detailed and more precise p2| . Moreover, we shall consider a large number of fragmentation processes: p — > tt^, 
7t^~ — > tt^ 1 , K + — > 7r ± , and ir + — > If , and the results of our calculations can all be compared with existing data. To 
our knowledge such data have never been used before to confront model calculations. That the valon-recombination 
model is able to reproduce those data, as we shall show in this paper, should therefore be regarded as meeting a 
higher demand than fitting the charge asymmetry of ir + p collisions discussed in [jj) . 
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Another reason for revisiting the soft processes in hadronic collisions is that the recent development in heavy-ion 
collisions at high energies presents some urgency to understand better the basic processes of particle production at 
a more fundamental level. To understand the formation of dense matter, it is necessary to understand first baryon 
stopping and pionization in pA collisions. We have recently done a comprehensive treatment of the problem of 
momentum degradation in pA collisions in the framework of the valon model . The emphasis is on the nuclear 
effects on the projectile. The parameters used for the valon and parton distribution functions have not been put to 
test on the hadronic collisions. Since those distributions have been updated even more recently ||, it is necessary to 
focus on the elementary processes, not only of pp, but also of meson-proton collisions, in order to check the valon- 
recombination model (VRM) to a degree never attempted before. It is toward that end that we apply the VRM to 
the soft processes in this paper. 



II. VALON AND PARTON DISTRIBUTION FUNCTIONS 



The valon model describes the hadron structure relevant for multiparticle production. The momentum distributions 
of the valons can be determined from the parton distributions at low Q 2 , as posted by CTEQ4LQ fll4| |. That was 
done very recently |9), resulting in excellent fits of the u and d quark distributions. The parameters describing the 
valon distributions turn out to be quite different from the old values based on imprecise muon and neutrino data of 
the 70's ||, although the formalism of the model remains the same. Furthermore, the assumption of the symmetric 
sea used previously has been lifted so that the parton distributions in the valons are now also very different. Here we 
give a summary of the valon and parton distributions functions, the details of which can be found in j^] . 

In a proton the 3-valon distribution is 

GuuD(yi,y2,V3) = g P (ymTvi S(yi + 2/2 + 2/3 - 1), (1) 

where yi is the momentum fraction of the ith valon (y being never used for rapidity in this paper), and 

g p = [B(a + 1, p + l)B{a + 1, a + /3 + 2)] _1 , (2) 
B(m, n) being the beta function. The single valon distributions are obtained by integration 

G v {y)= [ dy 2 f dy 3 G UUD (y,y 2 ,y 3 )=g p B(a + l,p+l)y a (l-y) a+ P+\ (3) 



G D (y) = J d Vl J dy 2 GuuD(yi,y2,y) = g P B{a + l,a + l)y\l-y) 2a+l . (4) 

The new values of a and (3 are found to be || 

a = 1.76 and (3 = 1.05, (5) 

which are significantly different from those used in Refs. [l3| due to various theoretical assumptions and limited 
experimental data. 

Since each valon contains one and only one valence quark of its own flavor, the valence quark distributions are 
convolutions of the relevant valon distributions with the valence (non-singlet) quark distributions in the valons. In 
terms of moments we have simple products 

u v {n) = 2G u {n)K NS {n) 1 (6) 
d v (n) = G D (n)K NS (n), (7) 



where 

Gu, D (n)= I dyy n -'Gu, D (y), 1*) 



l 

n — 1/ 







and similarly for q v (n) in terms of q v (x) where q stands for either u or d. The value of Q 2 is set at l(GeV/c) 2 for 
application of the model to low-p-r processes and will not be exhibited explicitly. From Eqs. (||) and (Q) we have 

Gu(n) =B(a + n,a + /3+ 2)/B(a + 1, a + + 2) (9) 
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G D (n) = B(fi + n, 2a + 2)/B(j3 + 1, 2a + 2). (10) 
The n dependence of KNs( n ) is determined in || and fitted by the following parameterization 

I 3 A 

KNs(ri) — exp — Cjw' , u = ln(n — 1), (11) 

V j=° / 

where c 3 = 0.753, 0.401, 0.0962, and 0.0555, for j = 0, 1, 2, 3, respectively. 

For the sea quark distributions in the valons we have found that the SU(2) symmetry has to be broken to fit the 
CTEQ4 data, and that the favored quark distributions (u in U and d in D) are suppressed relative to the unfavored 
quark distributions (u in D and d in U), consistent with Pauli blocking. They are denoted, respectively, by Lf(z) and 
L u (z), which are different from the s-quark distribution L s (z) and the gluon distribution L g (z). The distributions of 
these partons in the proton are given, in terms of the moments, simply by sums of products 

ti = 2GuL f + G D L u , (12) 
d = G D L f + 2GuL u , (13) 



8=[2G U + G D )L„ (14) 

g = (2Gu + Go) L g (15) 
where the dependences on n have been omitted. The functions Li(n) are parameterized as follows fl 
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In Li{n) = -^bfvP, u = ln(n - 1) (16) 

3=0 

(i) 

where b) are given in Table I. 



TABLE I: Coefficients in Eq. @ 
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4.12 


2.2 


0.2 


0.18 


u 


3.07 


1.5 


0.08 


0.05 


s 


4.21 


1.6 


0.1 


0.02 


9 


0.98 


1.0 


0.05 






For meson-initiated reactions we need the valon distributions in mesons, which can be determined from the quark 
distributions in mesons on the basis of the universality of the quark distributions in valons. From experimental 
data on Drell-Yan and prompt photon production in n^N collisions, the parton distributions in the pion have been 
determined in Ref . , using the parameterization 

xq v (x) = A v x a ' {1 - xf (17) 

for the valence quarks with a 1 = 0.64 ± 0.03 and fj 1 — 1.11 ± 0.04. In the valon model for the pion that distribution 
is related to the valon distribution by 

xq v (x) = / dyi / dy 2 G K (yi, V2) K NS (x/yi) , (18) 
Jo Jo 

where Kpjs( z ) is identified with the same in the proton problem. Since the two valons in pion are symmetrical, we 
use the valon distribution 



G w (2/1,2/2) = c/tt (2/12/2) 7 $ (yi + 2/2 - 1) 



(19) 
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where = l/Bi^y + 1,7 + 1)- We convert Eq. ( p"8| ) to the moment form 

q v {n)=G*{n)K NS (n), (20) 



where G*(n) involves only one parameter 7, and ifjvs(n) is the same as in Eq. ( |ll| ) on the basis of the universality 
of the valon structure, independent of the host hadron. The corresponding moment of Eq. (|l7| ) can trivially be 
calculated. Both are shown in Fig. 1, where a good agreement between the two is achieved by the choice 

7 = 0. (21) 

This value is also consistent with the one used in Ref. j^] . Thus g v = 1, and the single valon distribution in a pion is 
simply 

G%y) = 1. (22) 

The flat distribution of the valons in pion is a result of the fact that the constituent quarks are much more massive 
than the pion so they are tightly bound, resulting in large uncertainty in the valon momentum fraction. 

The situation with the K meson is somewhat different, since the kaon mass is higher, and the constituent quarks 
have unequal masses. The valon distribution has a form similar to that of the pion in Eq. (fL9h 



G K (2/1, 2/2) =g K y a 1 yl6{y 1 + y 2 -l) (23) 
where qk = 1/ B(a +1,6+ 1). The average momentum fractions of the two valons are 

V1.2 = J dy 1 dy 2 yi.2G K (2/1,2/2) (24) 



so their ratio is 



yjy 2 = (a+l)/(b + l). (25) 



Since the average velocities of the two valons are the same as that of the host kaon, their momenta should be 
proportional to their masses, i. e., 2/1/2/2 ~ ra\j jras- Taking the constituent quark masses of u- and s-types to be in 
the ratio 2:3, we get from Eq. ( p5| ) 

6=(3a+l)/2. (26) 

Thus we are left with one unknown parameter in Eq. (^3|). Since the parton distributions in kaon are not known, wc 
cannot determine that parameter as in the pion case. It will, however, be determined by the soft production data of 
kaon-initiated collisions later in Sec. 5. 

Knowing the valon distributions in mesons results in our knowledge of the recombination functions for the formation 
of the same mesons. Since we define the recombination function in the invariant phase space, we have from ( |l9|) and 
(§) 

(x llX2 ,x) = ^5 (Tl + El-i) (27) 
Rk {x x ,x 2 ,x) = g K Q+1 6+1 S p + ^ - l) (28) 



with a and b being constrained by (|26J). In this paper we do not consider nucleon production, so R p is not needed 
(see pi). 



III. PROTON FRAGMENTATION 



In a pp collision the soft production process in the fragmentation region is treated in the VRM as one in which 
the proton bag is broken by the collision and the valons become clusters of partons. The central idea of the model 
is that the probability for detecting a meson at large x is higher for a q and a q at lower Xi to recombine than for a 
q or a diquark at high x' to fragment. The distributions of q{x\) and q(x 2 ) depend on the hadron structure. Under 
the assumption that the collision process does not significantly perturb the parton distributions outside the central 
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interaction region (i.e., Xi > 0.1), the valon model can provide a sensible link between hadron structure and parton 
distributions. Since glueballs have never been seen, the gluons hadronize by first converting to qq pairs, thereby 
enhancing the sea. Downstream the quarks and antiquarks dress themselves and become the valons to be recombined 
in forming the produced particles. Since the dressing process does not change the total momentum of a quark, the 
meson momentum x is simply the sum of x\ and X2 of the q and q, and the probability of hadronization is determined 
by an integration of the enhanced q{x\) and q(x 2 ) weighted by the recombination function. This is an s-channel 
description of the fragmentation process, the initial formulation of which is given in Ref. §. What we have now are 
the new distributions || and the more extensive low-pr data |ll| not available to be considered in || . 
Quantitatively, the invariant distribution for pion production in the proton fragmentation region is 



dx 



H„(x) = [ — — (xi,x 2 )Rk (xi,x 2 ,x) 
J Xi x 2 



(29) 



where F v (x\, x 2 ) is the invariant distribution for q at x\ and q at x 2 , and is a convolution of the valon distribution in 
a proton with the quark distributions in valons, whose sea quarks are enhanced. In Fig. 2 we show the subprocesses 
giving u and d quarks by schematic diagrams, which represent 



F„+ (x 1 ,x 2 ) 



dy\dy 2 [G uv (y 1 , y 2 ) K (x 1 /y 1 ) L' u {x 2 /y 2 ) 

+2G UD (y u y 2 ) K {x 1 /y 1 ) L' u (x 2 /y 2 ) 
+ 2G UD ( yi ,y 2 ) L' u (xi/i/i) L' u {x 2 /y 2 )} 

12 



J dy 2G v (y) jjf fe) L' u 



where 



K(z) = K NS (z) + L' f (z). 



y-xi 

X-2 



y-xi 



12. 



(30) 



(31) 



L'f(z) and L' u (z) are favored and unfavored sea quark distributions in a valon, enhanced by gluon conversion, and will 
be discussed below. In Eq. ( |30] ) the symbol {. . .}i2 implies symmetrization in x\ and x 2 . Equation ( (30| ) can appear 
simpler in moment form, for which we define 



/■I pl—xi 

F 7r +(n ll n 2 )= / dxi / dx 2 x" 
Jo Jo 



ni —2n2 —2 



X T F„+(x l ,x 2 ), 



(32) 



so that (pOf) becomes 



F Tr +(n 1 ,n 2 )= G uu (n 1 ,n 2 ) + 2G UD (ni,n 2 ) K (m) L' u (n 2 ) 

+2G UD (m,n 2 )£;(ni)i; (n 2 ) 
+2Gu {n x + n 2 - 1) [k (m,n 2 ) L' u (na)} 

+G D (ni + n 2 - 1) $yL' u (ni,n 2 ) L' u (n 2 )| 

where 

K(n u n 2 )= [ dzz ni - 2 {l- z) n2 - l K(z) 
Jo 

and similarly for L' u (ni,n 2 ). For n~ we have 

F 7T -(n 1 ,n 2 )= Guu (ni,n 2 ) +2G UD (ni,n 2 ) L' u (nx) L' u (n 2 ) 

+2G UD (m, n 2 ) (m) ^ (n 2 ) 
+2G;7 (ni +n 2 - 1) (rii,n 2 ) ( n 2)| i2 

+G D (m + n 2 - 1) {i^ (m,n 2 ) l4 (n 2 )} . 



(33) 



(34) 



(35) 
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To use these in Eq. (E9f) we note that the moments of H'^(x) — x 3 H n (x) is, using Eq. (12 



1 r l fl-Xl 

n— 2 rr/ , 



H4(n)= / dxx^H'^x)^ / dxi / dx 2 ^(xi,x 2 )(xi+x 2 ) Tl 
J o Jo Jo 



(36) 



where the sum over n\ and n 2 is restricted by the constraint n± + n 2 = n. The double moments of the 2-valon 
distributions are 

Guu (ni,n 2 ) = / dy 1 dy 2 yi 1 ~ 1 y 2 l2 ~ 1 GuuD (2/1, 2/2,2/3) 
-1 ri-i/i 



/3 



3p / «2/i/ "2/22/1 2/2 (1-2/1-2/2) 
Jo Jo 

= g p B (ni + a, n 2 + a + (3 + 1) B (n 2 + a, [3 + 1) , 



(37) 



Gao (ni,n 2 ) = / dy-idysvT 1 y% 2 1 G UUD (2/1,2/2,2/3) 
= .g p £? (ni + a, n 2 + a + /3 + 1) B (n 2 + /3, a + 1) , 



(38) 



where g p is given by Eq. (g). 

For the enhanced sea quark distributions, L'f(z) and L' u (z), we recall that the quiescent sea as probed by electroweak 
interaction is described by Lf(z) and L u (z), given in Ref. For either a U or D valon, the momentum fractions of 
all its partons add up to one, so we have a constraint on the n = 2 moments 



i^ s (2) + 2 L f (2) + L u (2) +2L S {2) + L g (2) = l 



(39) 



The various terms above correspond to (for U valon, say) u valence, uu sea, dd sea, ss sea and gluons. For hadronization 
we consider the enhanced sea where the gluons are completely converted into the uu and dd sectors, a scenario which 
we refer to as the saturated sea. Note that we do not let the gluons be converted to the ss sector due to the higher s- 
quark mass. Such a restriction should be relaxed in the case of nuclear collisions because of substantial Pauli blocking 
in the non-strange sectors. Thus we write 



From these two equations follows 



^jv S (2)+2 L' f {2) + L' u (2) +2L S (2) = 1 



L g {2) = 2 L'{2) + L' u {2) - L f {2) - L u {2) 



(40) 



(41) 



Assuming that the enhancement factor f q is the same for favored and unfavored quarks, i.e., L'j u (2) = / g L/ jtl (2), we 
get 



f q = l 



L g {2) 



L f (2) + L u (2) 



(42) 



From Table I, the values of In 1^(2) are found to be —4.12, —3.07 and —0.98, for i = f, u, and g, respectively. One 
thus obtains 



/, = 3.99. 

We assume that this enhancement factor applies uniformly at all z so that we have 

L 'f,u( z ) = f q LfA z ), or L' f ,u(n) = /<?£/,«(«)• 



(43) 



(44) 



This assumption is made mainly for the sake of simplicity and to avoid introducing undetermined parameters. It 
cannot be expected to be valid at very low z where the density is high, but that is outside the region of applicability 
of VRM anyway. 
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The two last quantities in Eqs. ( |33| ) and ( p5| ) that remain to be specified are K(n\,n 2 ) and L' u (n\,n2). They 
both are defined by integrals of the typed given in (|34|) in terms of Kns(z), L'Az) and L' u (z), the latter two being 
proportional to Lf(z) and L u (z). Since the single moments Kjys( n )j Lf(n) and L u (n) are given in Eqs. (|ll| ) and 
(p"6|), we consider a direct approach to calculating K(n\,ti2) and L(, (711,712). Denoting Kns(z), Lf{z) and L u (z) 
collectively by J(z), we determine the double moments 



J(m,n 2 ) = / dzz ni - 2 {\ - z]^- 1 J(z) (45) 
Jo 

by expanding (1 — z)™ 2_1 in powers of z up to O (z™ 2 ), since J(z) is small as z — » 1, and the binomial coefficient c m 
decreases rapidly with to. Thus to a very good approximation we can write 

J(ni,n2)= £(-l) ro ( K("i+™)> (46) 

m=0 ^ ' 

from which we can calculate K(ni,n 2 ) and L' u (ni,n 2 ). 

Having specified all the terms in Eqs. ( |33| ) and (|35"|), we can now calculate -F 7r + (^1,^2), which in turn are used in 
( |36| ) to determine H^±{n). The inversion to H'^ ± (x) [and then trivially to H 7T ±(x)] involves a method that has been 
discussed in Ref. |l3| and is summarized in the Appendix. We show our results in Fig. 3. The solid lines are from 
our calculations with normalizations adjusted to fit the data jl2|, since the data are for the dimensionful Ed 3 a/dp 3 
at fixed pt while our model is for the dimensionless xdN/dx, integrated over all pr- The shapes of the calculated 
inclusive distributions involve no adjustable parameters. Evidently, they agree very well with the x dependences of 
the data. We regard this result as evidence in support of the VRM. 

IV. PION FRAGMENTATION: tt+ -> vr 

In this section we consider the non-diffractive inclusive distribution of 7r + — > n~ in the fragmentation region of the 
pion beam. By charge conjugate invariance it should be the same as for n~ — > 7r + , if there is no contamination by 
target fragmentation. Since there are no valence quarks of the projectile that end up in the produced pion, one naively 
would expect dN/dx to be strongly suppressed at high x. Yet experimental data (T^j indicate that its x dependence 
is rather similar to that of p — > 7r + , in which a u quark is shared between p and n + . Thus simple fragmentation 
of a valence quark or diquark in a string model cannot account for this similarity. In the VRM because the valon 
distribution in the pion is different from the ones in the proton, with the former giving larger average momentum 
fraction of the valons than the latter, large- a; 7r~ in 7r + can arise from the higher-momentum valons. In this section 
we show this behavior quantitatively. 

For 7r + — > 7r~ the invariant distribution H n +_> % - (x) has the same formal expression as Eq. (^9|) except that for 
p fragmentation is to be replaced by F ir +_^ T -. In Fig. 4 we show the schematic diagrams for UD valons going to du 
quarks. The moments for F ir +_ tn - can then be written down by inspection 

(num) = Gud (ni,n 2 ) L' f (ni)L' f (n 2 ) + i^(ni)L^(n 2 ) 

+ [G£ (m) + Gl (n 2 )] {^(ni,n 2 )i u (n 2 )} i2 . (47) 

We note that u is favored in U, and d is favored in D. Because we assume that the structure of valons is universal in 
all host hadrons, the enhanced sea of a valon in the pion is the same as that of a valon in the proton; thus L'j u is the 
same as in the previous section. 

Since G^g (2/1,2/2) = 5 (yi + 2/2 — 1) due to 7 = in Eq. (|l9|), we have 

Gl D {n 1 ,n 2 )=B{n u n 2 ), (48) 
GId = !M- (49) 



The calculation of F V +^ IT - (ni, n 2 ) is therefore straightforward. Using it in Eq. (36) yields i? lr +_ >7r - (n), and by 
inversion we get 11^+^^ - (x). 

The result is shown in Fig. 5 by the solid line. The dotted line is the fit of the data given by the experimental 
paper [12], where the x > 0.6 points are excluded in order to avoid contamination from resonance decay products. 
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Our result agrees very well with the experimental parameterization of the data by ~ (1 — x) m where m = 3.37 ± 0.09 
for px = 0.3 GeV/c at beam momentum 100 GeV/c. 

The lack of K in Eq. (|47|) compared to <J33[) and (]3q), corresponding to no valence quarks, is compensated by the 



fact that G^jq and in (|4|) and (ffl) are not as damped at high m as G^, G V VD and G v v D are in (|37|), (|38|), (||) 



and ([L0|), respectively, for a and /3 being as big as in (|5j). That is, in n + we have harder valons and softer d and u, 
while in p we have softer valons and harder u and d. 

The normalization of the calculated distribution in Fig. 5 is again adjusted to fit, since the data are for inclusive 
cross section at fixed px- Nevertheless, the agreement of the x dependence with data is the second piece of support 
one can infer for VRM. 



KAON FRAGMENTATION: 



In Sec. 2 we have discussed the form of the valon distributions in a kaon with one undetermined parameter to be 
fixed. Now we consider kaon initiated reactions and determine that parameter by fitting the inclusive distribution of 

K+ -> 7T ± . 

From Eq. we have 



G* s (ni,n 2 ) = B (m + a, n 2 + b) / B{a +1,6+ 1), 



(50) 



G£ (m) = B (m + a, b + 1) /B{a + 1,6+1), 



(51) 



Gf (n 2 ) = B (a + 1, n 2 + 6) /B(a + 1,6+1), 



where a and 6 are constrained by (|26|). In Fig. 6 we show the schematic diagrams for (a) 
It then follows that 



J J W+-*+(ni,"a) = G$ § (m,n 2 ) \K( ni )L^{n 2 ) + L' u { ni )L^{n 2 ) 

+G* (n 1 ){k(n 1 ,n 2 )L' u (n 2 )} i 



(52) 

7r + and (b) — > 7r _ . 



+G 



A' 



(najjL^m.naJL^na)} 



.12 



(53) 



where L"(n 2 ) is the moment of the enhanced, unfavored quark distribution in the S valon, and will be discussed 
below. Similarly, we have for K + — ► 7r~ 

F K+ ^- (ni,n 2 ) = G^ 5 (m, n 2 ) f^(m)L^(n 2 ) + L' u { ni )Ll{n 2 ) 



+G« (n 1 ){ZUn 1 ,n 2 )L' / (n 2 )} i 
+Gf (n2){^(ni,na)Z2(n2)} 12 



(54) 



For gluon conversion in the 5 valon we again consider the saturation of the uu and dd sectors of the sea, but none 
in the ss sector because of the higher s-quark mass. Since the non-strange sectors are both unfavored, the momentum 
carried by the sea is not the same as that in a non-strange valon. Thus we write before gluon conversion 



K NS {2) + 4Lf (2) + 2Lf (2) + If (2) = 1 



(55) 



where Lf(i = u,s,g) refer to the S valon. Equation ( p5[ ) differs from ( |39| ) for a D valon only in that the favored 
sector dd is replaced by the unfavored sector dd. It means that there is a redistribution of the momenta in the gluons 
and sea quarks, which, we assume, takes the simple form 



Lf(2)=cZ l (2), 



where c is a constant. After gluon conversion Eq. (55) becomes 

X JVS (2)+4^(2) + 2Lf(2) 



(56) 
(57) 
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If we define L^(2) = f s L s s (2), we obtain 



f s = 1 + L s g {2)/AL S U {2) = 1 + L S (2)/4L„(2) 



Using as before L u (2) = e~ dm and L g (2) = e - u ys , we get 



(58) 



(59) 



f s = 3.02. 

To determine £"(2) in terms of L u (2) we need, in addition to f s , the value of c in Eq. (|5^). It follows from Eqs. (|3S|), 
@ and (||) that 



c = 1 



2 L u (2)-L f (2) 



(60) 



4L U (2) + 2L S (2)+L g (2) 
Since Z s (2) = e -4 21 from Table I, we get 

c = 0.9 (61) 
and finally L" a {2) = f s cL u (2) = 2.72L U (2). Again, extending this proportionality to all n 2 , we have 

Z£(n a ) = 2.72L u {n 2 ). (62) 



We now have all the quantities in Eqs. (52) and fl54|) to calculate Fk+^„± (m , n 2 



There is one adjustable parameter in our calculation. It is a in Eq. (23), b being constrained by (26). After 
Fk +— >7r± ("l; n 2) are determined, we use a formula similar to Eq. (pq), i.e., 



H' K +^± 0) = Y\ -^—,F K+ ^± (ni + 2, n 2 + 2), 



(63) 



to calculate H' K+ ^ ±(n). Then by inversion to H' K+ ^ 7r± (x) as before, we can finally obtain H K +^ 7r ±(x) — 
x~ 3 H' K+ ^ n± (x). We adjust a to fit the data, which are shown in Fig. 7 for pt = 0.3 GeV/c. Again, because 
the data are for fixed p?, we cannot predict the normalization, which is adjusted to fit. The best values of a and b are 



a = 1.0 



6 = 2.0. 



(64) 



The solid lines in Fig. 7 show the calculated results; the dotted lines are the experimental fits using the form (1 — x) m . 
Although the fits are not perfect, they are acceptable in view of the large error bars in the data. What is noteworthy 
is that we used only one normalization factor for both n + and 7r - production, and the calculated curves agree with 
the two sets of data in their relative normalizations. Thus the VRM has captured the essence of production in 
K + initiated reactions. 

Since a and b are now known, we can exhibit the valon distributions in a kaon. In Fig. 8 we show (y) and 
Gg(y). Note how the strange valon has larger momentum fraction than the nonstrange valon, on the average. It is 
because of the harder S valon, though softer non-strange quarks in it, that gives rise to the distributions H K +_ t7T ± (x) 



that are harder (decreasing more slowly with x) than either n + — > n or even p 



which has a shared valence 



quark. This is a remarkable affirmation of the importance of the structure of the hadron in the determination of the 
inclusive distributions of its fragments. 



VI. PION FRAGMENTATION: ty+ -> 



Since the valon distributions in the kaon have been determined in the previous section, we now know the re- 
combination function for the formation of kaon. Thus the calculation of pion fragmentation into should be 
st r aight forward . 

In Fig. 9 we show the diagrams for tt + — > K^. They are essentially the same as the ones in Fig. 6 for K + — > 7r ± , 
except for the replacement of S valon in the initial K + by D valon in the initial 7r + , and the replacement of the 
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appropriate non-strange o and q quarks in ir by the s and s quarks in K . The corresponding equations are 



therefore similar to Eqs. (p2|) and 

F 7T +^ K +{ni,n 2 ) = Gu B (n ll n 2 ) K{ni)L s (n 2 ) + L s {ni)L' u (n 2 ) 

+GI (m) {k(n u n 2 )L s (n 2 )} ^ 

+Gl(n 2 ){L' u ( ni ,n 2 )L s (n 2 )} u (65) 



Ftt+^k- (ni,n 2 ) = GIj B (m, n 2 ) Z^(ni)L s (n 2 ) + L s (n 2 )l4(ni) 

(ni){L^(ni,n2)L a (n2)} 
+G£ (n 2 ){ZUni,n 2 )Z s (n 2 )} i2 



(66) 



Note that L s (rii) is not enhanced because, as before, gluon conversion saturates the non-strange sectors of the sea. 
The other moments are as given in Sec. 4. 

Since the recombination function Kk is now different from R v , the inclusive distribution for production is 



= H K (x) = [ -^-—^-F^K (2:1,3:2) Rk (xi,x 2 ,x) 

J Xi x 2 

dx\dx 2 x\x\Fn-,K x 2 ) S (xi + x 2 — x) 



x dN K 
dx 

= 9KX 

have been used. If we define 

H' K (x) =x 6 H K (x), 

then the binomial expansion of (x\ + x 2 ) n leads us to the moments 



where Eqs. 



and 



H' K (n) = / dxx n - 2 H' K (x) = g K V 

JO r„ 1 



ni\n 2 \ 



F„-,k (ni + 3,n 2 + 4) 



(67) 



(68) 



(69) 



similar, but not identical, to Eq. (|36]). In view of Eqs. ( p5|) and (g6|), H' K ± (n) can now be calculated without any free 
parameters. 

Using the same procedure as before to obtain H K ± (x), we can determine the inclusive distributions for n + — > K . 
In Fig. 10 we show our results compared to the data [O at 100 GeV/c and px = 0.3 GeV/c. Again, the normalization 
is adjusted to fit, but only one normalization factor for both curves. The agreement between theory and experiment 
is very good, considering that no free parameter is used except for the overall normalization. The excess for x > 0.6 
can be attributed to the decay of K* whose K product would have higher x. The relative normalization of K + to 
K ~ is well reproduced by VRM. 



VII. CONCLUSION 

We have considered the production of mesons in the fragmentation regions of incident proton and mesons in \ow-px 
collisions. The data on p — » n^, tt + — * it~ , K + — * 7r ± , and 7r + — ► K ± have all been satisfactorily reproduced by 
the calculated results in the VRM. There is essentially only one free parameter that is connected with the valon 
distribution in the kaon. All other parameters specifying the structure of the proton and pion have been determined 
independently by fitting the data on hard processes. It is self-evident that by successfully reproducing the inclusive 
distributions of all the above reactions we have met the test of charge asymmetry in ir + p collisions. 

Apart from the details of the VRM and of the data in the fragmentation region, the fundamental themes that this 
work has affirmed are that the hadron structures are important and that hadronization proceeds through recombina- 
tion. The valon model effectively describes the hadron structure and the valon distributions provide the probability 
functions for recombination. It is conceivable that two valons may be regarded as a diquark whose fragmentation 
yields the produced mesons. Since the two valons are spatially distributed objects with non-vanishing relative mo- 
mentum, to describe their fragmentation by a fragmentation function adapted from jet physics seems hard to justify. 
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The VRM, on the other hand, is related more closely to the parton model and provides a natural s-channel description 
of the fragmentation process (from the point of view of the hadron) in terms of recombination (from the point of view 
of the quarks and antiquarks). 

We have not considered the production of baryons in this paper. However, the non-diffractive production of nucleons 
in pA collisions has already been treated in the VRM . The production of strange particles is worthy of further 
attention. Our consideration of 7r+ — > K is only a beginning, which shows that in hadronic collisions gluons do not 
convert to ss as effectively as to uu and dd. That situation must change in going from hadronic to nuclear collisions 
due to Pauli blocking in the non-strange sector. It seems that the VRM provides a natural framework in which to 
investigate the transition from strangeness suppression to strangeness enhancement. 
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APPENDIX A: APPENDIX 



We summarize here the method of inversion from the moments to the x-distribution function, originally proposed 
in Ref. Instead of making the inverse Mellin transform, which involves a complex contour, we exploit the 

orthogonality of the Legendre polynomials. First, we shift the variable to the interval < x < 1, and define 



= P e (2x - 1) 



so that 



dxg t (x)g m (x) 



I 



2t + l 



If we expand the distribution H'(x) in terms of gt{x) 

oo 

H\x)=Y i {2l+l)h l g t {x), 

then the inverse is 

he = / dxH' (x)ge(x). 
Jo 

These hp can be expressed in terms of the moments H'(n) if we express gi{x) as a power series in x 



fji: 



(x) =J2 a e 



(Al) 



(A2) 



(A3) 



(A4) 



(A5) 



where a\ are known from the properties of Pi(z). Thus from Eq. (A4) we have 



h/>=^ra}H'(i + 2), 



(A6) 



where H'(n) is defined in Eq. (|36|). It is now clear that our theoretical results in H'(n) can be transformed to H'(x) 
through Eqs. ( |A3j ) an d (|A6| ) once we have the coefficients a\. Furthermore, if H'in) becomes unimportant for n > N, 
then the sum in Eq. (A3) can terminate at N. 



To determine al, we make use of the recursion formula 



+ l)Pe+i(z) = {21 + \)zP t [z) - m-i(z) 



(A7) 



to infer through Eqs. (|Alj) and (A.5) 



4 



; [(2£-l)ati + (e-l)al 



2h 



(A8) 



4 - - - + V ~ l)4- 2 ], (A9) 

where I > 2, and 1 < i < I. For I < 2, we have 

og = l, a? = -l, a} = 2. (A10) 

With these we can generate all a^, so /i^ can be directly computed. The use of ( [A3] ) then yields H'(x). 

We have found that this method can give very accurate result in inverting H'(n) to H'(x) for N roughly between 
8 to 10, depending on how rapidly H'in) decreases with n. 
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FIG. 1: Moments of the valence quark distribution in a pion, q v (n), where the parameter 7 is chosen to fit the moments of 
distribution determined from the experimental data [p.5|. 
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FIG. 2: Schematic diagrams in the VRM for the production of tt and n in the proton fragmentation region; only the ud 
and du states are shown. 



14 





15 



10 



-a 



b 



1 



10 



Data at 1 OOGeV 
p T =0.3 GeV/c 
Calculated 




0.2 0.4 v 0.6 0.8 



FIG. 5: Inclusive distributions of tt + —> n . Data are from Ref. [|12|. The dotted line is the experimental fit; the solid line is 
the theoretical result in the VRM. 
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FIG. 6: Schematic diagrams in the VRM for (a) K + —* ud, and (b) K + — > du. 
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FIG. 7: Inclusive distributions of K+ -> ty+ and if+ -> vr". The data are from Ref. |l2| P L = 100 GeV/c and p T = 0.3 
GeV/c. The dotted lines are the experimental fits; the solid lines are the theoretical results in the VRM. 
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FIG. 8: The momentum-fraction distributions of the valons in a kaon: (a) U valon in solid lines, and (b) S valon in dotted 
line. 
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FIG. 9: Schematic diagrams in the VRM for (a) tt + — > us, and (b) n + — » sit. 
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FIG. 10: Inclusive distributions of tt+ -> K + , and (b) tt+ -> AT - . The data are from Ref. @ Pl = 100 GeV/c and p T = 0.3 
GeV/c. The dotted lines are the experimental fits; the solid lines are the theoretical results in the VRM. 



